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Abstract 

The Skyrme crystal, a solution of the Skyrme model, is the lowest energy-per- 
charge configuration of skyrmions seen so far. Our numerical investigations show 
that, as the period in various space directions is changed, one obtains various other 
configurations, such as a double square wall, and parallel vortex-like solutions. We 
also show that there is a sudden "phase transition" between a Skyrme crystal and 
the charge 4 skyrmion with cubic symmetry as the period is gradually increased in 
all three space directions. 
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1 Introduction 



The Skyrme model was originally proposed by Tony Skyrme in 1961 [lj. It is a theory of 
nuclear matter in which the fundamental building blocks are pion fields. It was later shown 
[21 13] that it can be regarded as a low-energy approximation to QCD, which becomes more 
accurate in the large A^-limit, where N c is the number of quark colours. The (topologically 
non-trivial) solutions of this theory are called skyrmions. There is a special way of arranging 
them in order to obtain the smallest known value of the energy per baryon number seen 
so far [3], known as the Skyrme crystal. This is an infinite, triply-periodic, arrangement 
of half-skyrmions j5j EH [HI IHJ HUJ. In order to study Skyrme crystals, we impose periodic 
boundary conditions, with period L, on the Skyrme field in all three directions. The 
skyrmions are therefore defined on a 3-torus T 3 . 

Skyrme crystals were originally proposed by Klebanov [5j as a model for dense nuclear 
matter, such as that found in neutron stars. At the time, the behaviour of two well- 
separated skyrmions was already known [lj - an important feature being that these are 
maximally-attracted when one is rotated with respect to the other by 180° about a line 
perpendicular to the line connecting them. Klebanov believed that an interesting extension 
of this idea would be to have an array of skyrmions, where any skyrmion would be attracted 
by its nearest neighbours. A graceful way of achieving this, Klebanov showed, would be to 
arrange the skyrmions in a simple-cubic lattice with appropriate rotations applied to all 
the nearest neighbours of a chosen skyrmion. Moreover, if the skyrme fields are to have the 
correct periodicity in all three directions, they must have symmetry elements that combine 
both spatial transformations as well as isospin transformations acting on the individual 
fields (more on this later in the introduction). 

Klebanov also showed that there is a minimum in the energy per baryon number of the 
lattice for a certain value of the period L. In other words, there is a preferred size of the 
fundamental cell. Later, Goldhaber and Manton showed J6] that there is a phase transition 
from a low-density simple-cubic lattice of skyrmions to a high-density body-centred lattice 
of half-skyrmions. However, it has since been shown by Kugler and Shtrikman [9j ED] and 
by Castillejo et al. [8] that the lowest energy per baryon configuration is that of a (high- 
density) half-skyrmion phase corresponding to an initial (low-density) face-centred cubic 
(fee) array of skyrmions. 

An fee array is one in which skyrmions with standard orientation are placed on the 
vertices of a cube and more skyrmions are placed on the face centres, but this time rotated 
by 180° about an axis perpendicular to the face. Such a configuration produces 12 nearest 
neighbours (to a particular skyrmion), which are all in the attractive channel. If the origin 
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is fixed at the centre of one of the unrotated skyrmions and the skyrme fields are given by 
^(x) = (<I>i(x), $ 2 (x), $ 3 (x), $ 4 (x)), then the fee configuration would have spatial and 
isospin symmetries, as was alluded to above. The generators for these symmetries are listed 
in Table 1 0j: 



Table 1: Symmetry generators of an fee array of skyrmions 



Transformation Name 


Spatial Transformation 


Isospin Transformation 




(x u x 2 ,x 3 ) (-x u x 2l x 3 ) 


($1, $ 2 , $4) $2, $3, $4) 


K- 




($1, $ 2 , $3, $4) ($2, $3, $1, $4) 




(x u x 2 ,x 3 ) -> (x u x 3 , -x 2 ) 


($1, $3, $4) ~> ($1, $3, ~$2, $4) 




(x u x 2 , x 3 ) (x 1 + L/2, x 2 + L/2, x 3 ) 


($1, $3, $4) (-^1, ~$2, $3, $4) 



The transformations listed here are 1Z\'. a reflection in the x\— axis, 1Z\ _: a "negative" 
three- fold rotation about the diagonal that goes from the origin to the opposite corner of the 
cube (defined as "1"), 1Z\ _\ a "negative" four-fold rotation about the X\ — axis, and 7i+, 2 +: 
a positive L/2— translation in both the x\ — axis and the x 2 — axis. Note that these are a 
subset of the possible transformations that can be carried out on an fee lattice. For example, 
one can also have a "positive" three- fold rotation along the same diagonal, TZ\ + , given 
by the transformation: (xi,x 2 ,x 3 ) (x 3 ,£i,x 2 ), ($1, $ 2 , $ 3 , $ 4 ) ($ 3 , $1, <£> 2 , $ 4 ). 
Another possible transformation could also be a "positive" four- fold rotation about the 
x 3 -axis, Kj +: (x u x 2 , x 3 ) (-x 2 , x u x s ), ($1, $ 2 , $ 3 , $4) (-$2, $1, $3, $4). 

There is an additional symmetry unique to the high-density phase of half-sky rmions. 
Its generator is given by: 



Table 2: Additional symmetry generator of the high-density half-skyrmion phase 



Transformation Name 


Spatial Transformation 


Isospin Transformation 


T l+ 


(x u x 2 , x 3 ) -> (xi + L/2, x 2 , x 3 ) 


($1, $ 2 , $3, $4) (-^i, $ 2 , $3, -^4) 



Note that this transformation involves a chiral SO (4), rather than just isospin SO (3) 
rotations displayed in Table 1 above and it can replace the 7i+, 2 + transformation since 
that can be achieved through successive applications of 7i+ and 75+. 

In this letter, we will see that the field behaves in an interesting way as one changes 
the period of the crystal in different ways, for different directions. For example, if we start 
with the Skyrme crystal and then increase the period along all three space dimensions 
in the same way, one gets the familiar picture of the cubically-symmetric charge Q = 4 
skyrmion. We also show that the energy density of the Q = 4 skyrmion displays certain 
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similarities previously seen in the context of Skyrme chains, provided we have large L x ^ y 
(where L x = L y are the periods in the x\ — and x 2 — directions) and small L z values. Skyrme 
chains are solutions of the Skyrme model, which are periodic in one space dimension. It 
has been shown pj] that soliton chains generally have constituents in the form of vortex- 
ant ivort ex pairs. These emerge when the period is small compared to the natural soliton 
size - a feature that we verify when the period in the z-direction is small compared to the 
other two space dimensions. When the period increases, the constituents tend to clump 
together, a feature that is also verified here. 

A double Skyrme sheet [12], is also seen to emerge at small L x ^ y and large (and effectively 
infinite) L z values. It takes the form of a square lattice, an object analogous to the 
hexagonal "Skyrme domain wall" solution [13]. However, our system is periodic in all 
three directions, which means that the vacuum value on both sides of the Skyrme sheets 
(±oc in the z-direction) is unique. 

Finally, we describe what happens as one increases the period simultaneously in all 
three directions, starting with the Skyrme crystal, and show that there is a rapid tran- 
sition between the Skyrme crystal and the Q = 4 skyrmion with cubic symmetry. We 
show evidence which suggests that this is a second-order phase transition with an order 
parameter given by the period L x ^ y ^ z of the configuration (where L x = L y = L z ). 

2 The Skyrme Crystal 
2.1 Background 

The static energy density of the Skyrme model is given by: 

S = -^Tr(L^) - -^Tr([Lj, Lj][L i} Lj]) , (1) 

where = U^dU / dx\ x l = (x, y, z) are the spatial coordinates, and the field U (x l ) is an 
SU(2)-valued scalar field. The pions are grouped into this scalar field as follows: 

U = $ 4 + i$i<n , (2) 

where $^ is a triplet of pion fields, Oi is the triplet of Pauli matrices, and $4 is an additional 
scalar field determined through the constraint: WU = QpQp = 1. 
The energy E is defined as 

E = — ^— / Edxdydz . (3) 
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The scalar field U is a map from ]R 3 , compactified at infinity, to 5 3 , the group manifold 
of SU(2). Skyrme identified the degree of this map, a topological invariant, with the baryon 
number, which is given by 

Q= Qdxdydz, (4) 



where 

Q = ^- 2 e l3k Tx{UL 3 L k ) (5) 

is the topological charge density. The energy ^ satisfies the Faddeev-Bogomolny lower 
bound [H] E > Q. 

An analytic approximation for the fields of the Skyrme crystal was proposed in j8j. It 
takes into account the SO (4) chiral symmetry as well as the SO (3) isospin symmetries. 
The fields are expressed as follows: 

$ 4 = cic 2 c 3 , (6) 

$1 = —Si ( 1 — — H — ^ j and cyclic permutations, (7) 

where Si = sin(27rx z /L) and q = cos(27rx 2 /L). It is a good approximation to the actual 
minimal-energy solution. 



2.2 Changing the periods L x = L y and L z 

In what follows, we take periodic boundary conditions in all three space directions - the 
periods will be specified in the relevant sections. The lattice spacings in the x, y, and z 
directions are given by h X) h y) h z and the number of lattice points are given by n X) n y) n z , 
yielding side-lengths L XiViZ = h x ^^ z *n x ^^ z . We use a first-order finite-difference scheme and 
implement a full 3-dimensional numerical minimization of the energy using the conjugate 
gradient method (cf. [15]). 

Note that there is a numerical error associated with the finite lattice spacing. The 
way we approximate the errors in the energy, for a given configuration, is by comparing 
its topological charge (using numerical methods) with the "true" value of its charge. We 
assume the same errors in energy and charge since similar finite-difference methods are 
employed in calculating each of these values. We have noticed that higher differences in 
the lattice spacing in each space direction yield higher errors. 

The initial condition that we start with for minimization is the approximate skyrme 
crystal of eight half-skyrmions, namely ^ and 0. After being minimized, the period in 
each direction for this initial configuration is then changed in a certain way (described in 
the relevant section) and then re-minimized. 



5 



2.2.1 From a 4-skyrmion to a square 2- wall 

The first case is the one for which the period in all three directions is large for the initial 
configuration and then the L XjV periods are reduced gradually. The initial period is given 
by L x ^ y ^ z = 7.05 (see Fig. [T]) and then reduce L x = L y by 1 each time. 

We start by increasing the periods L x = L y = L z for the minimized Skyrme crystal 
from L Xj y jZ = 4.7 to 7.05 by increasing the value of n x<g , z , for a certain value of h Xj y jZ , and 
thus producing Fig. 1(a). Afterwards, n Xi y jZ is kept constant and L Xj y jZ is reduced by 1 each 
time. 

One can see that the translation symmetries, 77+/_, (where i = {1,2,3} and in either 
direction +/— ) of the high-density Skyrme crystal are broken in Fig. 1(a), which is the 
(5 = 4 skyrmion, whereas all reflection symmetries, 7t^, three-fold rotations 72^ 2 +/-j an< ^ 
four- fold rotations lZf + /_ remain unbroken - characteristic of cubic symmetry. In Fig. 
l(b)-(f), 7i+/- and 72+/- symmetries are regained (possibly due to the fact that we are 
"squeezing" the configuration in these directions), whereas 75+/- remains broken only for 
Fig. 1(b) and Fig. 1(c), due to the fact that these figures are not extended throughout 
the whole period in the z— direction. The reflection symmetries IZi are unbroken and the 
three-fold rotations 72^ 2 +/- are broken in Fig. l(b)-(f). The four-fold rotations TZ^ + ^_ 
remain unbroken throughout, whereas TZy 2 +/- are broken in Fig. l(b)-(f). 

Note that a double Skyrme sheet configuration [12] emerges in Fig. 1(d). The separation 
between the sheets, which is calculated by measuring the distance between the energy 
density peaks as a function of z, remains constant as L X y is reduced further and L z is kept 
constant at L z = 7.05. However, the energy decreases from a value of E — 4.36 ± 0.09 at 
L x , y ,z = 7.05 down to a minimum of E = 4.20 ± 0.05 at L Xjy = 4.05 and, finally, increases 
to E = 5.25 ± 0.10 at L XtV = 2.05. 

The preferred configuration for this Skyrme sheet is to have L x = L y . For instance, if 
we fix L y = 4.05 and vary L x away from this value in either direction, we notice that the 
energy increases. 

2.2.2 From square 2- walls to vortices 

We now start with a large L z period, which will be reduced, and keep L x<y = 4 constant. 
The starting point is L z = 7, which is then reduced by 1 each time. In this section, 
n x = n y = n z are kept constant, producing different h z values each time L z is changed. 

The initial condition (eqs. ^ and Q) is first minimized for the periods L XiV = 4 and 
L z = 7, producing the double Skyrme sheet configuration seen in Fig. 2(a). The only 
broken symmetries associated with this configuration, which is seen to persist in Fig. 2(b) 
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(a) L x y =7.05, L z =7.05 (b) y =6.05, L.=7.05 (c) l_ x y =5.05, L z =7.05 




Figure 1: Energy density isosurfaces with surface value given by 0.5 * £ max , where £ max is 
the maximum value of the energy density (all isosurfaces in subsequent figures have the 
same surface value). The first isosurface corresponds to the Q = 4 skyrmion for the periods 
L x , y ,z = 7.05. Each successive picture has L XiV reduced by 1 and L z is kept constant at 
7.05. 

and Fig. 2(c), are 1Zy 2 +/- an d 7^i/ 2 +/-- These are regained in Fig. 2(d), which is a 
(non-minimal) Skyrme crystal configuration. As the L z period is decreased further, we 
notice the appearance of vortex-like structures, which are periodic in the z— direction (Fig. 
2(e) and Fig. 2(f)). The symmetries associated with these configurations are the same as 
those in Figs. 2(a)-(c) and, in fact, the translation symmetry in the z— direction, 75+/-, 
becomes continuous for Fig. 2(f), since there is no longer any noticeable z— dependence. 
This ties into the subject of Skyrme chains, which has been explored in [llj. 

The energy of the isosurfaces decreases from a value of E — 4.22 ±0.05 at L z = 7, down 
to a minimum of E — 4.17 ± 0.04 at L z = 5 (where the isosurfaces are still in the form of 
a double square wall), then back up to E = 4.89 ± 0.05 at L z = 2. 

2.2.3 From square 2-walls to vortices, via crystal 

In this section, we show how one can transform a double Skyrme sheet configuration into 
the 4- vortex configuration discussed in the previous section, by changing both L x = L y 
and L z (rather than just L z ) ) and going through an intermediate, minimal, Skyrme crystal 
state, as can be seen in Fig. [3j The initial period is L x ^ y = 2.7 and L z = 10.7. The former 
is then increased by 1 and the latter decreased by 3, two consecutive times, producing Figs. 
3(b)-(c). The changes in the periods are then swapped, increasing L x ^ y by 3 and decreasing 
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(a) L xy =4, L z =7 (b) L xy =4, L z =6 (c) L <y =4, L z =5 




Figure 2: Energy density isosurfaces corresponding to the Q = 4 double Skyrme sheet 
configuration, where L x ^ y = 4 remains fixed and where L z = 7 is reduced by 1 in each 
successive picture, transforming the double skyrme sheet into a (non-minimal) Skyrme 
crystal and thereafter into a parallel 4-vortex structure. 

L z by 1, for two consecutive times, producing Figs. 3(d)- (e). Here, n x = n y = n z are kept 
constant, producing different h x = h y and h z values, each time the periods are changed. 

We start by minimizing the approximate Skyrme crystal (eqs. ^ and ([7])) for the 
periods L x ^ y = 2.7 and L z = 10.7, producing the double Skyrme sheets discussed in the 
previous sections (with the same symmetries) in Fig. 3(a). As the periods are changed, as 
described above, the Skyrme sheets still persist in Fig. 3(b), and change into the (minimal- 
energy) Skyrme crystal in Fig. 3(c) (with all its associated symmetries as described in the 
introduction). As L x ^ y increase and L z decreases further, the Skyrme crystal changes into 
a 4-vortex structure, Figs. 3(d)-(e), which have less of a z— dependence than the ones in 
Figs. 2(e)-(f), respectively. 

The Skyrme crystal has the minimum energy, with E — 4.13 ± 0.04, followed by the 
square 2-walls at L XiV = 3.7 and L z = 7.7 with E = 4.26 ± 0.06. The final picture, 
the 4-vortex configuration with the smallest L z — value, has the highest energy, with E = 
4.77 ± 0.17, which follows from the fact that Skyrme chains "prefer" to be closer together 
in the x, y— directions [IT] . 

Thus far, we have changed the periods in the three space directions in different ways - 
by decreasing L x ^ y (or keeping it constant) and keeping L z constant (or decreasing it) or 
changing all three periods at different rates at the same time. We now turn to the case 
where we increase (or decrease) all three periods simultaneously at the same rate. 
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(a) =2.7. L z =1 0.7 (b) L, y =3.7, L z =7.7 (c) l_ x y =4.7, L z =4.7 




Figure 3: Energy density isosurfaces, which show a transition from a pair of Skyrme sheets 
to a 4-vortex configuration. Subsequent pictures have L XjV and L z changing in such a way 
that their values are swapped halfway through the transition. 

2.2.4 From the Skyrme crystal to the Q = 4 skyrmion 

It turns out that an interesting feature of the Q = 4 system is uncovered when one starts 
from the minimal-energy Skyrme crystal configuration at L x ^ y ^ z = 4.7 and then increases 
the periods simultaneously until the individual half-skyrmions coalesce. 

These half-skyrmions clump in a sudden fashion when one perturbs the system in a 
certain way. We perturbed it by removing one lattice site through the middle of the 
configuration in all three directions, thereby "squeezing" the half-skyrmions together in all 
directions. Therefore, in order to see them coalescing, one needs to increase the periods 
very gradually. This appears to happen in the range L x ^ y ^ z = 6.07 — 6.09 as can be seen in 
Fig. [3J where the first picture has periods L XiViZ = 4.7 and the following three correspond 
to the range just mentioned^] Compare this with Fig. [5| where the periods are increased 
in larger steps. 

The sudden merging of the half-skyrmions can be visualized in a different fashion by tak- 
ing the maximum value of the difference in the energy densities of the fields £[Qp{xi)] under 
a certain symmetry transformation, in this case a translation of L x /2 in the x— direction, 
which is half the size of the fundamental cell, and dividing by the maximum value of the 
energy density, i.e. A x = (£[$p(x 3 )] - <?[$/3(x j/ )]) max /^ max , where x 3 ' = x J - L x /2. We 

x It is worth mentioning that the value of this range changes when the bin size is changed - we used the 
value n x ^y^ z = 32. However, the general feature that the half-skyrmions coalesce rapidly remains the same. 
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(a) L xy z =4.7 (b) L x y z =6.07 (c) ^ y z =6.08 (d) ^ y z =6.09 




Figure 4: The energy density isosurface of the Q = 4 Skyrme crystal at = 4.7. The 

periods are then increased gradually from L XiViZ = 6.07 to L XiViZ = 6.09. 



(a) L x y z =4.7 (b) L xyz =5.88 (c) L >; =7.05 




Figure 5: The first figure corresponds to the energy density isosurface of the Skyrme crystal 
for the periods L XiViZ = 4.7, the second figure has periods L XjViZ = 5.88, and the third one 
has L Xj y jZ = 7.05 

then plot this as a function of the period L XjViZ - the reason we do this is that, as soon as 
the half-skyrmions start to coalesce, they will no longer be L x /2— periodic. For the Skyrme 
crystal, this difference is seen to be essentially zero and it then starts to increase as the 
half-skyrmions begin to coalesce as can be seen in Fig. [6] 

This jump in the asymmetries of the crystal is analogous to a phase transition in 
thermodynamics. There is a sudden transition from a crystalline phase, which has more 
symmetries, such as chiral SO (4) symmetries for the Skyrme crystal, to a phase with less 
symmetries, such as SO (3) isospin symmetries for the Q = 4 cubic-shaped skyrmion. We 
expect the transition to become more pronounced as the number of bins increases - as can 
be seen in Fig. 6, when n Xj y iZ is increased from 32 to 36, which starts to resemble a step 
function. It should be noted that the Ai values are independent of which L XiViZ value one 
starts with, which hints to a lack of hysteresis in the system. Extending the analogy with 
phase transitions, this lack of hysteresis strongly suggests a second-order phase transition 
(cf. [IE]), with no latent heat, and with the period L x ^ y ^ z as the order parameter. 
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Figure 6: Difference in the value of the energy density of the Skyrme fields under an 
L x /2— translation in the x— direction divided by £ max , for different values of the period 
L x , y ,z- Note that as the bin size decreases (or as n increases), the plot begins to look more 
like a step function. 

3 Concluding Remarks 

We have seen that changing the periods of the Skyrme crystal, away from the values which 
minimize its energy, has allowed us to find new configurations - such as the double square 
lattice of "Skyrme sheets" - as well as unearthing unforeseen connections with Skyrme 
chains. The Skyrme sheets are produced when the period L x ^ y is small compared with L z . 
Our results [12] have shown that these are more energetically favourable than the double 
hexagonal case, which makes them a prime candidate for use as building blocks for large- 
charge configurations. If one then goes to the other limit of large L Xjy and small L z one gets 
a series of four vortex-like charge Q — 1 objects, which is related to the subject of Skyrme 
chains [llj. As the periods are changed and we go through these different configurations, 
we have seen that these generally merge and separate in the directions in which one is 
changing the period, for example, the merging of the half-skyrmions in the z direction as 
the period is decreased along this direction. 

We have also noticed that the phase transition from the Skyrme crystal to a charge 
Q = 4 skyrmion is likely to be second-order, due to the lack of hysteresis in the system, 
and with an order parameter given by the period L x ^ y ^ z . 

Note that the deformations considered here have all been along one or more of the edges 
of the Skyrme crystal, keeping it rectangular. An interesting extension, which we plan to 
investigate, would be to see if any new structures emerge when one applies more general 
affine transformations, such as a stretch along a diagonal of the crystal. 
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